Current driven rotating kink mode in a plasma column with a non-line-tied free end 
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First experimental measurements are presented for the kink instability in a linear plasma column 
which is insulated from an axial boundary by finite sheath resistivity. Instability threshold below the 
classical Kruskal-Shafranov threshold, axially asymmetric mode structure and rotation are observed. 
These are accurately reproduced by a recent kink theory, which includes axial plasma flow and one 
end of the plasma column that is free to move due to a non-line-tied boundary condition. 
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The current driven kink instability is a magnetohydro- 
dynaniic (MHD) instability which affects current carry- 
ing plasmas in Nature and laboratory. The kink mode 
structure and stability condition are strongly dependent 
on the system geometry and the boundary conditions 
(BCs). KruskalQ and Shafranov (hereafter referred 
to as KS) considered first the ideal MHD stability of a 
cylindrical plasma column with magnetic field compo- 
nents (0, Bg, Bz) using cylindrical coordinates (r,9,z). 
For an infinitely long (equivalent to periodic axial BCs) 
column, they obtained a linearly unstable helical kink 



mode of structure ^ = e 



i{d+2Trz/L) 



when the plasma cur- 



rent In exceeds the Kruskal-Shafranov limit 



IKS 



{2^Ya'B,/{^ioL) 



(1) 



where a and L are, respectively, the radius and length 
of the current channel, and ^ is the displacement of the 
plasma column from the equilibrium position. The KS 
theory has been quite successful in predicting the behav- 
ior of toroidal plasmas for which the periodic BCs yield a 
proper accounting for the finite length of the system. In 
linear systems, however, substantial deviations from KS 
predictions can result from different axial BCs. 

The importance of the BCs has long been recognized^, 
0] and is of particular relevance to the stability of line- 
tied flux ropes in space physics (c./. Ref. and survey 
Ref. 0), and astrophysical jetsQ- In recent years, there 
has been a renewed interest in the stability of a line-tied 
plasma column in laboratory devices (see Refs. 0, 0| 
and references therein). The kink stability of a plasma 
column with line-tied ends has been investigated in linear 
devices, where line-tying is attributed to the presence of 
highly conducting end plates[ll]|. and in open systems to 
a local discontinuity for the Alfven velocity that forms a 
virtual boundary around the svstem|l2l 

In this Letter, we experimentally investigate the exter- 
nal kink instability in conditions where one end of the 
plasma column is line-tied to the plasma source, and the 
other end is not line-tied and therefore free to slide over 
the surface of the end-plate. The latter BC is a result 
of plasma sheath resistance that insulates, at least par- 
tially, the plasma from the end-plate. Compared to the 
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FIG. 1: Experimental setup with diagnostics and coordinate 
system. (1) Vacuum vessel; (2) plasma gun; magnetic probes 
in (3) an axial array and (4) an azimuthal array; (5) external 
coils; (6) external anode. On the right, an axial cut near 
the external anode. The fast camera (7) is located at the 
midplane and views the plasma column along the x direction. 
The bi-dimensional magnetic probe (8) measures {5Bx, SBy) 
at the edge. Also schematically shown is the plasma column 
whose end rotates at the external anode. 



line-tied case, we find significant differences in the kink 
mode structure and lower critical current for the onset of 
the kink. The axial flow velocity, the direction of the ax- 
ial magnetic field, and the rotation frequency of the kink 
mode are all correlated. The experimental results agree 
with the predictions from a recent theory of the external 
kink instability for a slender plasma column [l^. 

The experiments are conducted in the Reconnec- 
tion Scaling Experiment (RSX) which uses electrostatic 
plasma guns to generate the plasma^^. Figure ^ shows 
the experimental setup with a view of the plasma gun 
located at z = and radially inserted into the center of 
the RSX cylindrical vacuum vessel (4 m length and ra- 
dius Tjuaii — 0.2 m). A hydrogen arc plasma is ejected to 
form a cylindrical plasma column in a constant, uniform, 
axial magnetic field of B^ = 0.012 T generated by exter- 
nal coils. Fig. CKS). Electron density and temperature 
have Gaussian profiles with half-maximum radius ~ 2 
cm such that rg r^aii- Central electron density and 
temperature are rieo = 1 — 3 x 10^^ m~'^ and Teo = 5—14 
eV, and axially decrease away from the gun. The axial 
velocity Vz « 3 — 5 x 10^ m"-'^ of the plasma has been 
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FIG. 2: Time histories of (a) the plasma current; (b) az- 
imuthal magnetic field Bg at (z = 0.48 m, 6 = n/2, r = 0.15 
m). Expanded view of the temporal evolution of 5 Be from 
azimuthal (c) and axial (d) arrays of magnetic probes during 
the stationary phase of the discharge. 



estimated by solving the momentum balance equations 
for dominantly axial flow constrained by experimentally 
measured axial gradients in pressure and density. A cur- 
rent is driven in the plasma by negatively biasing the gun 
with respect to an external anode (0.05 m^ stainless steel 
plate) which is electrically isolated from the vacuum ves- 
sel. The external anode can be positioned at different 
axial locations z = 0.3 —3m allowing the length L of 
the current-carrying plasma column to be varied. 

The MHD activity is monitored by two arrays of mag- 
netic probes that include a total of seven coils inserted 
in the vacuum vessel at Vp = 0.15 m to measure the az- 
imuthal magnetic field Bg. In the axial array, Fig. ^-(3), 
four magnetic probes are positioned at 6* = 7r/2 (top of 
the vessel) and axially located at z = 0.14, 0.48, 0.62 and 
0.76 m. In the azimuthal array. Fig. 0(4), the four mag- 
netic probes are located at z = 0.48 m and azimuthally 
equispaced by 7r/2. One magnetic probe is shared by the 
two arrays. Signals are digitized at 20 MHz. 

In Fig. Efa), the time evolution of the plasma current 
measured at the external anode is shown for a discharge 
with L = 0.92 m. The plasma current starts at t = 
when the bias is turned on (arc discharge starts at t = — 1 
ms) and after sa 150 //s reaches a stationary phase which 
lasts for w 1.4 ms. The total plasma current Ip ~ 320 
A during the stationary phase is determined by the bias 
voltage and can be varied in the range Ip = 0.05 — 1 kA. 
Figure |2tb) shows the time evolution of Bg as measured 
by the magnetic probe located at z = 0.48 m, 6 = tt/2. 

For an expanded time window during the stationary 
phase, measurements of the perturbed azimuthal mag- 
netic field SBg from the azimuthal. Fig. OJc), and the 
axial. Fig. |2Id), arrays of probes show oscillations with 
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FIG. 3: (a) Amplitude of the perturbed azimuthal magnetic 
field, SBe, and (b) rotation frequency of the mode during the 
stationary phase as a function of Ip. (c) Axial structure of SBe 
close to the instability threshold for Ip/Icrit ~ 1-2 showing 
non line-tied BC at the external anode. The dot-dashed line 
shows the structure of SBg as predicted from theory. 



TO = 1 azimuthal periodicity at the frequency / ~ 50 
kHz. These oscillations are observed when Ip exceeds 
a current threshold Icrit (see below), but are not de- 
tectable for Ip < Icrit- The observed instability growth 
time tg « 4 ^s is of the order of the axial Alfven time. 

The m = 1 azimuthal structure is consistent with a 
plasma column rotating as a rigid body at frequency /. 
The phase of the signals from the axial array increases 
linearly along the z direction. Fig. El^d), indicating a 
rotating helical current channel. The kinked deformation 
is a right (left) handed helix if • i?^ > (< 0), as 
expected from a paramagnetic kink^^l- The direction of 
rotation reverses when reversing the external magnetic 
field Bz- In both cases, the mode rotation is such that 
the helix always screws into the external anode. Both 
the rotation frequency and the average amplitude of the 
SBg oscillations at the dominant frequency during the 
stationary phase scale linearly with Ip as shown in Fig. 
|3{a,b) where the plasma current is varied in the range 86 
A < /p < 190 A for a plasma column length L = 0.92 m. 
The plasma current corresponding to the limit 5Bg 
corresponds to the kink stability threshold which for 
these data is Icrit = 70 ± 7 A. The rotation frequency at 
the stability threshold is fcrit 28 ± 3.5 kHz, Fig. ^h). 

Figure l^c) shows the axial structure of the to = 1 
mode close to criticality for Ip w 1.2 x Icrit- While a 
line-tied BC holds at the gun, the large amplitude of 
the mode at z > 0.76 m can only occur if the plasma 
slides over the external anode and is therefore not line- 
tied. This is confirmed by images of the plasma column 
at the external anode as shown in Fig. ^ Ha emission is 
imaged with a Cooke DiCam Pro intensified cameraflij 
that views the external anode edge perpendicularly to 
the z axis, Fig. ^-(7). In Figs, ^a-b), two fast-gated 
(200 ns gate ^ Alfven time) images show respectively 
the plasma column in the equilibrium position and at 
the maximum displacement during the rotation. 

The observed sliding of the plasma column over the 
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FIG. 4; Optical and magnetic measurements near the exter- 
nal anode demonstrate PNLT BC. As the plasma gun shoots 
from right to left, a fast gated image of visible emission shows 
(a) the equilibrium, yeq, and (b) the displaced, ymax, posi- 
tion of the plasma during the rotation, (c) Measured optical 
displacement (l = \ymax — yeq\ for different Ip values, (d) 
Magnetic field perturbation SB = {SB^ +SByy^^ versus IpCl 
compared with predictions from Eq. @ (dashed line). and 
SBx,y are measured 0.01 m in front of the external anode. 



external anode surface is not consistent with a line-tied 
BC, which might be presumed at the external anode 
for time scales less then a 100 — 200 /iS resistive dif- 
fusion time. Deviations from line-tying may appear if 
the plasma is insulated from the external anode by a 
finite plasma sheath resistance. The degree of insula- 
tion is measured by the ratio k of the Alfven transit to 
inductive decay times of the current We estimate 
K w 15 from the formula k w {cs /vz)[c/ (atOpi)]'^ pl^'^ 
where Cg = {Te/rriiY^'^ , (3^ = rieTf./B'^ and uipi is the 
plasma frequency (derivation in Ref. foi' typical 

plasma parameters a « 2 cm, Tg « 10 eV, rig « 10^^ 
m^'^, Vz/cs ~ 1. For the present experiments, we con- 
clude that K 3> 1 and the plasma is insulated from the 
conducting boundary by finite sheath resistance. Ryu- 
tov and coauthors [lj| have shown that for k ^ 1 the 
BC can be formulated as {dC/dz + i{ko/2)()\z=L = 
where fco = Bg/^aBz) = /io/p/(27ra^i3z) with Cartesian 
complex displacement C, = (,x{z,t) + i(,y{z,t). In the fol- 
lowing, we will refer to this condition as to perfect non- 
line-tying (PNLT) BC and we will show that it holds at 
the external anode, after using an equivalent condition 
that can be directly compared with experimental data. 

With a thin plasma column approximation (a <C 
fwaii J L) and assuming that r = is the equilibrium 
position, the Cartesian components of the perturbed 



magnetic field at r 



^ a can be expressed as 



SB,,y = ia^/rl)B,i±dCx^y/dz+koCy^x)- Using the PNLT 
condition, in the vicinity of the external anode the per- 
turbed magnetic field can be expressed as 



SB = (SBl 



-5Bl)^'^ 
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We use magnetic measurements combined with images 
near the external anode to test Eq. (0), equivalent to 
the PNLT condition. Images such as those in Fig. 01a,b) 
determine the maximum displacement ymax and equilib- 
rium yeq positions of emissivity along a line parallel to the 
external anode surface and axially spaced 0.01 m in front 
of it. The image in Fig. Efb) is synchronized to a mag- 
netic probe at the same axial location but azimuthally 
separated by 7r/2 to provide the position of maximum 
column displacement. In Fig. EJc), the displacement 
Cl = IVmax — Veql is shown for a series of discharges with 
100 A < Ip < 500 A. In Fig. 2||d), the magnetic pertur- 
bations {SBx, SBy) are measured at rp ~ 0.17 m and 1 
cm in front of the external anode using a bi-dimensional 
magnetic probe. In Fig.^d), the experimental data show 
excellent agreement with the PNLT predicted scaling in 
Eq. (0) (shown as a dashed line). We therefore conclude 
that the PNLT BC applies to the external anode. 

For a PNLT BC at z = L and with axial plasma 
flow velocity Vz, the eigensolution for the MHD equa- 
tion of motion (Eq. (57) in Ref. [H) is C = C(e*'=^^ - 
gjfc z-jgie-iwt ^i^Y^ axial wave numbers k± = ±(7r/2L)[l=F 
(1 - M2)1/2] where M = Vz/{V2va) is the Alfven Mach 
number, and va is calculated using the average density 
fie — (2/a^) J^ne{r)rdr. In Fig. ^c), the axial struc- 
ture of SBg from experimental data for Ip w 1.2 x Icrit 
(squares) and from the theoretical eigensolution (dot- 
dashed line) are compared. The good agreement suggests 
a robustness of the eigensolution for slightly super-critical 
regimes. The total axial phase shift at criticality (not 
shown here) is also consistent within the experimental 
uncertainties with the expected eigensolution for Mach 
numbers in the range M = 0.25 — 0.4. 

The critical current in this case is 

Icru = ii^^^&(l - Aiy/- ^ ^(1 - M^/- (3) 

showing that the kink mode is unstable at half the KS 
current for vanishing flow and smaller currents when flow 
exists. The eigenfrequency lu at the criticality has a finite 
real component 



i?e(c.) = -^^M(l-M^) 



2N1/2 



(4) 



where CL = [UL,t)'+^y{L,tf]^/\ 



The perturbed plasma column rotates at frequency 
Re{uj) ~ 27r/, constant along the column. The rota- 
tion is driven by axial plasma fiow from the gun to the 
external anode along the helically-perturbed column as 
well as a Doppler shifted k± ■ Vz, such that the helix al- 
ways screws into the external anode, consistently with 
the experimental observations. 

To investigate the applicability of Eqs. ® - Q, figure 
shows (a) the critical current and (b) the rotation fre- 
quency at the criticality as a function of the inverse col- 
umn length 1/L which is in agreement with a predicted 
linear scaling. 
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FIG. 5: For three different column lengths (L — 
0.48,0.92,1.38 m): (a) critical current Icrit (b) rotation fre- 
quency at the criticality fcrit as a function of the inverse col- 
umn length 1/L. The grey region indicates predictions from 
KS theory. Dotted lines show linear fits. 



The system comprising Eqs. ® -O) provides a theo- 
retical relationship between Icrit, fcrit, a,, M and the av- 
erage plasma density (or Alfven speed). This was solved 
using experimental measurements of the electron den- 
sity profile at criticality together with deduced scalings 
= 59/L [A] and fc„t 2.4 x W^/L [Hz] from the 
data in Fig. El Wc obtained a radius a = 2 ± 0.2 cm 
which is in agreement with the experiment (see Fig. Ej) 
and Mach number M = 0.25—0.4. The corresponding ax- 
ial flow velocity « 3 — 4.8 x lO'' m~^ is also consistent 
with our independent estimate. In Fig. [SJ predictions 
from KS theory using the computed a are compared to 
the experimental data. 

For the external kink, the radius a is expected at the 
position where conductivity, magnetic diffusion time, and 
axial current density Jz{r) are negligibly small. This is 
tested using profile measurements at the criticality. Fig- 
ure ini shows a Bcnnet profile least-squares fit of Jz (r) 
at the instability threshold, as obtained from azimuthal 
magnetic field measurements at z = 0.48 m for L = 0.92 
m. The previously determined radius a ~ 2±0.2 cm is lo- 
cated at the edge of the profile where Jr(a) « 0.05x Jz(0). 
Similar conclusions follow from the magnetic diffusivity, 
derived from the measured electron temperature profile 
and the observed instability growth time tq. 

It should be emphasized here that the measurements 
in this Letter only characterize the global structure of 
the external kink mode. Although some of the obser- 
vations (e.g., axial structure and rotation frequency of 
the mode) show some resemblance with measurements of 
internal drift-Alfven waves in finite-/? plasmas with ax- 
ial current [it), additional internal measurements will be 
needed to investigate the presence of these waves in our 
experiment. 

In summary, we have presented first experimental re- 
sults for the current driven kink instability in a plasma 
column with one free end that is not line-tied to its 
boundary due to a finite sheath resistance. Mode struc- 
ture and instability threshold arc significantly different 
from predictions from the Kruskal-Shafranov theory and 
are accurately reproduced by a recently developed kink 




FIG. 6: Radial profile of current density (solid line) from 
measured azimuthal magnetic field (squares) at the criticality. 
The vertical grey region indicates the error bounded estimate 
of a at the edge of the current channel. 



theory. The finite rotation of the mode, which is also ob- 
served in other laboratory dcvices|lll Il3|. is due to the 
plasma flow along the helically kinked plasma column. 

Though not addressed here, the effect of a non line-tied 
end should be important for the relaxing kinked plasma 
column. A destabilizing Jg x Bz force coexists with a 
stabilizing curvature force due to axial field line bending. 
The free end allows shifts with reduced curvature and 
restoring force possibly resulting in a helical saturated 
state of larger displacement than the line-tied case. Nu- 
merical simulations are being implemented for the non- 
linear evolution of the kink mode. 
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